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Abstract 

Recent developments in the synthesis of colloidal particles allow for control over shape and inter¬ 
particle interaction. One example, among others, is the so-called “Mickey Mouse” (MM) particle 
for which the self-assembly properties have been previously studied yielding a stable cluster phase 
together with elongated, tube-like structures. Here, we investigate under which conditions a fluid of 
Mickey Mouse particles can yield phase separation and how the self-assembly behaviour affects the 
gas-liquid coexistence. We vary the distance between the repulsive and the attractive lobes (bond 
length), and the interaction range, and follow the evolution of the gas-liquid (GL) coexistence 
curve. We find that upon increasing the bond length distance the binodal line shifts to lower 
temperatures, and that the interaction range controls the transition between phase separation and 
self-assembly of clusters. Upon further reduction of the interaction range and temperature, the 
clusters assume an increasingly ordered tube-like shape, ultimately matching the one previously 
reported in literature. These results are of interest when designing particle shape and particle- 
particle interaction for self-assembly processes. 
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I. INTRODUCTION 


One of the fundamental questions in condensed matter physics is to understand the 
relation between the microscopic and macroscopic properties of a system. The answer to 
this question is also closely related to the possibility of creating new routes to fabricate 
novel functional materials. A variety of fabrication protocols has been devised over the past 
years, but more recently the focus has been shifted towards encoding information for the 
self-assembly into the basic colloidal building blocks P-jS]. This paradigm is now known as 
self-assembly or “colloidal LEGO”. Beforehand, it is unknown how the encoded information 
will manifest itself at the end of the self-assembly process, and it is even unknown what kind 
of encoded microscopic information will yield the desired macroscopic structure. Answering 
these questions is the key to using the self-assembly process of colloidal particles as a new 
method to fabricate functional materials. 

“Hard” particle systems, which are purely entropic, give insights into the role of the 
particle shape in the self-assembly process. It has been shown that entropy alone can give 
rise to a large variety of close-packed crystal structures, as well as liquid and plastic crystal 
phases [MS]. On the other hand, having particles explicitly attract each other along specific 
directions can be of advantage to target particular open structures as, for instance, diamond 
crystals and kagome lattices[I3HIS]- This has brought a considerable amount of interest 
to what are now called patchy colloids, particles that interact with each other only via 
specific spots located at the particle’s surface [Tbll^ . Directional interactions can be given 
to colloids for example via attaching complementary DNA strands to the surface, as in Refs. 

luaEii. 

Another way of inducing directional interactions between building blocks is to combine 
depletion forces and surface roughness asvmmetrv [22l - [25] . In fact, if a building block is made 
up of rough and smooth lobes, the addition of depletants to the solution will induce specific 
attractions between the smooth lobes of different particles, while producing repulsions be¬ 
tween the rough-smooth lobes and the rough-rough lobes PH ES]- This approach has been 
used to self-assemble dumbbells j26] and, more recently, “Mickey Mouse” (MM) particles - 
trimers with one big smooth lobe (also referred to as the “head” of the MM particle) and 
two smaller rough lobes (the “ears ”) [27] . In particular, it has been shown that the MM 
particles self-assemble, at interaction strength of ~ 9 — IO^bT into tube-like structures [??]. 
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Yet, with respect to isotropically attractive spheres like, for example, Lennard-Jones 
particles or Hard-Sphere Square-Well (HSSW) particles, that exhibit a gal-liquid (GL) phase 
separation, it is so far unclear whether MM particles can also undergo an analogous gas- 
liquid phase separation (colloidal poor-colloidal rich) and, if this is the case, what is the 
interplay of the latter with the depicted self-assembly scenario. The aim of this work is to 
address this question with Monte Carlo (MC) simulations. To do so, we follow the evolution 
of the GL binodal line for a series of models connecting the isotropic HSSW spheres to the 
MM particles. The path proceeds firstly by increasing the rough smooth (ear-head) bond 
length until we reach the experimental particle geometry, and secondly by decreasing the 
interaction range until the experimental one is reached, as described in Ref. |27j . 

The paper is organised as follows: we introduce the model and the path connecting 
HSSW spheres and MM particles in Sec. |TI| where we also discuss details about the different 


computational techniques we have used. In Sec. |HI[ we present our results on the phase 
separation of MM particles and the transition to self-assembly. We also show how the 
particle geometry affects the structure of the liquid. We summarise and discuss these results 
in Sec. lIVl where we also provide an outlook for future research directions. 


II. MODEL AND METHODS 


We begin with the geometry of the MM particle as shown in Fig. [T} Each MM particle 
is represented by an aggregate of three spheres, two small spheres (“ears”), of diameter 
(Te, which represent rough lobes interacting with steric repulsions and a third bigger sphere 
(“head”), of diameter u/i, which corresponds to the smooth lobe and plays the role of an 
attractive site. The interaction Uj j between a pair of MM particles i and j depends on the 
positional and orientational degrees of freedom r*, r^, Slj, fij, respectively, which we drop 
here to lighten the notation, and consists of an attractive and a repulsive contribution, 


Uij = ulf + 


( 1 ) 


The attractive part ufh acts between the larger beads of particle i and j, and is given by 
a square-well (SW) interaction. 




J /3s for (Th < < Acr/, 

1^0 for > \ah 


( 2 ) 
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where Pe = e/k^T, e < 0 , represents the interaction strength compared to the thermal 
energy, A is the interaction range, and is the center to center distance from the head of 
MM particle i to the head of MM particle j. 

The repulsive part is determined by the hard-core interactions between the heads 
{h) and the ears ( 61 , 62 ) of two MM particles, and reads 


= XI -ry^l) (3) 

,62 

with the hard-sphere (HS) interaction. 




{ 00 if In,ry ,31 < cr „,,3 
0 otherwise 


(4) 


where = {aa -\- crp) 1^2. with a,P = h, 61 , 62 - 

Finally, we set the ear-head size ratio to q = aejcFh = 0.85 and the angle between the 
directions of the ears to O^e = 90°, in order to match the values of the experiments in 
Ref. [2Z1- However, the center-to-center distance between the ears and the head, denoted 
hereafter as the bond length l/oh^ is allowed to change, and it will be used as a parameter 
together with the interaction range. This particle and interaction model will be referred 
to as MMSW (Mickey Mouse Square-Well system). Note that in the limit I/ah = 0 our 
particles correspond to HSSW particles. 

We focus on the effect of varying the ear-head distance I/ah of the particle and the inter¬ 
action range A on the gas-liquid (GL) binodal line. We compute the GL binodal using the 
Successive Umbrella Sampling (SUS) technique pH] together with the Histogram Reweight¬ 
ing method [291-132]. Since the SUS method has been discussed elsewhere [28]. here we limit 
ourselves to only recalling the working scheme. The quantity of interest in the GL coexis¬ 
tence is P{N), the probability that the system will be in a state with N particles at hxed 
volume V, fixed temperature T, and fixed chemical potential /r. This probability is uni- 
modal above the critical temperature and assumes a typical bimodal shape, displaying gas 
and liquid peaks, for temperatures below the critical one. Furthermore, it can be shown 
that the bulk chemical potential at coexistence, at a given temperature, will be such that 
the area of the gas peak and the area of the liquid peak are equal. In fact, this is all we 
need to know to determine the coexisting chemical potential via the histogram reweighting 

technique [2nH32]- 
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FIG. 1. Model and parameters considered in this work. An experimental MM particle is obtained 
by a SW sphere by first changing the distance l/cih between the rough and the smooth beads (ear- 
head distance) and then by decreasing the interaction range A. The red beads (rough protrusions) 
have diameter (jg, while the blue bead (smooth protrusion) has diameter ah- 

The probability distribution function can be computed by splitting the entire N range 
into overlapping windows of fixed size and by performing Grand Canonical Monte Carlo 
(GCMC) simulations in each window. During the simulation one keeps track of how many 
times the system has a certain particle number N. which lies in between the lower and the 
upper limit of the window. The probability function can be reconstructed via the “stitching” 
procedure, which reads [28]: 

P(N) ^ ^ ftp) H,(N) 

P(0) ft(0) ft(l) ft(JV -1) ’ 

where we have implicitly assumed a window size of 1. From a computational point of view, 
Eq. shows us that this scheme is inherently parallel since all the ratios can be estimated 
from independent simulations. This is very convenient from a computer cluster perspective. 

Once P{N) has been calculated for a hxed temperature and (irrelevant) chemical po¬ 
tential, the coexistence chemical potential can be obtained by reweighting the distribution 
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until the area under the gas peak equals the area under the liquid peak. The reweighting is 
carried out via the following equation [SHI El]: 

lnP(A^|/5/ii) = lnP{N\/3iiQ) + - iio)N (6) 

In such a way, one can construct the GL coexistence envelope for particles with fixed geom¬ 
etry and fixed interaction range. In the following, we study how the GL coexistence curve 
is affected by the particle shape and particle-particle interaction range and where and how 
the regime of phase separation changes to self-assembly. 

Once we obtain the GL coexistence curve, we attempt to estimate the critical point via 
a least-squares fit to the equation: 

p±-p, = A\T- PI ±Ib\T- t/ (7) 

which stems from the law of densities and the law of the rectilinear diameter [331 135l - l37| . 
Here, ± stands for liquid/gas and (5 = 0.325 is the exponent of the 3D Ising universality 
class. Such a procedure yields only a rough estimation and we stress that an appropriate 
determination of the GL critical point should involve extensive use of the finite-size scaling 
technique [381 - llT] . 

To explore under which conditions of particle geometry, temperature and interaction 
range the interparticle attractions are more important than the repulsions, we additionally 
compute the second virial coefficient normalised to the one of a system of hard spheres 
with diameter Uft,, = B 2 /B^^. If attractions dominate, the value of the second virial 
coefficient wilt be negative, white positive otherwise. The temperature at which the second 
virial coefficient vanishes, so called Boyle temperature marks the crossing from 

one behaviour to another. The definition of the second virial coefficient involves the Mayer 
function, which takes into account the interaction between two MMSW particles [S^ ISH] 

fij = exp [-Puij] - 1 (8) 

The integral of the Mayer function over all possible positions and orientations of the two 
MM particles yields the second virial coefficient [i^ EH] 

B2(kBT/s) = -^ J dqidq2fi2 (9) 

where dq^ = dvi dfli represents integration over the particle’s positional and orientational 
degrees of freedom. Eq. can be computed via Monte Carlo integration. We place the 
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MMSW “1” with fixed orientation in the center of a box with volume V = (5ah)^ at a 
given temperature k-QT/e. We then generate a number W 0(10®) of random positions 
and orientations for the MMSW “2” in the same volume, and for each configuration k we 
compute the value of the Mayer function between particle “1” and “2”, /(^ 2 - Then, the 
second virial coefficient can be estimated asjSl 05] 

= = ( 10 ) 

fc =0 

where (/ 12 ) represents the average of Mayer function over all the configurations. We repeat 
this computation for different temperatures, particle shapes and interparticle range to locate 
regions where attractions prevail with respect to repulsions. 


III. RESULTS 

We have evaluated the GL coexistence curve for the MMSW model by varying the ear- 
head bond length I/ah from 0.1 to 0.57, corresponding to the experimental MM particles [27]. 
and the particle-particle interaction range A from 1.5 to 1.02. Following this path, we 
can connect the experimental MM system, for which the self-assembly behaviour has been 
investigated in Ref. 1271 , all the way to the HSSW limit, for which the GL coexistence has 
been intensively stndied [35l - [371 06] . 

A. Shift of the binodal line with the particle shape 

With a fixed interaction range A = 1.5, we change the particle shape by progressively 
increasing the bond length distance l/oh between the ears and the head of the MMSW 
particle. The highest value of bond length investigated in this work is set by the geometry 
of the experimental MM particle we want to match and reads l/uh = 0.57. We observe that 
the change in the MMSW particle shape affects the position of the GL coexistence curves, 
as seen in Fig. 

We quantify the change in the binodal curves, by computing the critical parameters as 
function of the bond length l/uh by means of a least square fit to Eq. Going from 
the HSSW to the MMSW, the critical temperature k^Tc/e shows a linear decrease upon 
increasing the bond length I/ah, as can be seen from Fig. In Fig. the critical critical 


7 



• = 0.10 ▼ I/Uf^ = 0.20 ■ = 0.30 * l/u^ = 0.40 ■ = 0.57 


FIG. 2. Gas-Liquid coexistence curves for MMSW particles. We provide both the temperature 
/ceF/e-density pa\, and temperature feeF/e-packing fraction 0 representation. The packing frac¬ 
tion is calculated as (/) = Nvp/V, with Vp the particle volume as discussed in Appendix The 
labels stand for different increasing values of the bond length 1. The full red dots indicate the 
location of the critical points by Eq. 

density pc is also seen to decrease monotonically with the bond length, while for the critical 
packing fraction (pc we observe a non-monotonic behavionr. 

Shifts of the critical temperature to lower values have been reported in different contexts. 
In fact, it is well known that attractive colloidal particles confined between two parallel plates 
display a shift of the critical point to lower temperatures because the confinement reduces 
the cohesive energy of the system[13-EI]. More recently, the concept of “limited valence” 
has been introduced in the context of patchy colloids [521 - I56] . The idea is that the valence of 
a particle (defined as the maximum possible number of bonded nearest neighbours) affects 
the location of a system’s critical point. Number of attractive patches |n2HSl|, patch surface 
cover age [nSl ES], patches with tunable attr actions [371 EH], ^re all examples of controlling a 
particle’s valence. Our model suggests that changing the distance between the attractive and 
repulsive spheres is another way of achieving such a control. In other words, we can think of 
the bond length as a “knob” for limiting the valence of the MMSW particle. Interestingly, 
similar trends for the critical temperature and critical density were recently reported for a 
system of attractive dimers interpolating between SW-SW and HS-SW in Ref. [3H] . 
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FIG. 3. Critical temperature k^Tc/s (full red dots) as a fuuctiou of boud leugth The greeu 

dashed hue shows the liuear fit to the data for the critical temperature /e = S + Al/a^^ 

with B = 1.37 aud A = —1.66, iudicatiug a liuear relatiou betweeu the boud leugth “kuob” aud 
the system’s critical temperature. 



FIG. 4. Critical deusity pc (full blue squares) aud critical packiug fractiou (pc (full cyau diamouds) 
as a fuuctiou of boud leugth I/ah- The solid curves are guides to the eye. 


The secoud virial coefficieut at the critical poiut has beeu suggested as possible measure 
of a particle’s valeuce iu Refs. [[55l [59]]. By usiug Eq. 10, we calculate the secoud virial 


coefficieut B 2 at the critical poiut [k^Tc/e for each value of the boud leugth I/ah- Our results 
are showu iu Fig. where we have uormalised the values by the secoud virial coefficieut of 
a hard sphere of diameter ah- Upou iucreasiug the boud leugth, the value of B 2 is seeu to 
decrease (become more uegative), which suggests that the limited valeuce picture could also 
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hold in the system at hand. 



FIG. 5. The reduced second virial coefficient at the GL critical point as a function of the ear-head 
bond length I/ah- The red line is a guide to the eye. 

It is interesting to investigate the structure of the MMSW coexisting liquid phase as 
function of bond length, by keeping the volume fraction hxed at (f) — 0.27. We do this by 
calculating the radial distribution function for all the MM attractive beads (heads) 
the results being shown in Fig. for bond length I/ah = 0.1 and l/ah — 0.57, respectively. 
We wish to remark here that we ignore the anisotropy of the particle and average over all 
possible particle orientations. Although both curves display features of the HSSW radial 
distribution functions, such as the presence of a hard-core diameter at r/ah = 1, and a 
discontinuity at r/ah = A = 1.5, we see that the peak positions for the case l/ah = 0.57 
hardly resembles the ones for l/ah = t).l. In fact, especially the hrst peak of the radial 
distribution function moves to higher values oir/ah, specihcally from r/cT/i = 1 for //ct/j = 0.1 
to = 1.5 in the case oil/ah = 0.57, and the same can be seen for the second and third 
peak. Clearly, even in the case of large interaction range A = 1.5, the presence of the ears 
of the MM particle has an influence on the fluid structure, pushing the heads further away 
from each other on average. 


B. Interaction range-driven transition to Self-Assembly 

After having analysed how the particle shape affects the GL phase separation, we now ask 
ourselves how this behaviour ultimately transforms into self-assembly, which is the behaviour 
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FIG. 6. Bulk head-head radial distribution functions g^^{r)ol the coexisting MMSW liquid at 
packing fraction (j) = 0.27 for bond length I/ah = 0.1 (magenta curve) and I/ah = 0.57 (green 
curve), as function of the distance r. The repulsive beads increase on average the spacing between 
the attractive beads. 

experimentally observed for the MM particles. In fact, the interaction range of the experi¬ 
mental MM particles system, i.e. A ~ 1.02, is much smaller than the one so far considered 
(A ~ 1.50). To this end, we perform a more extensive study of the second virial coefficient, 
where we map B 2 for several values of the interaction range, down to the experimentally 
relevant range, in the whole temperature k-oT/e and bond length plane I/ah- This gives 
us an overview of the conditions under which the attractions overcome the repulsions, and 
guides us in the choice of regions of the parameter space which are interesting to simulate. 

Fig. shows our results for the calculation of B 2 through the entire temperature U-qT/e 
bond length //(j/i plane for different values of the interaction range A = 1.50,1.40,1.30,1.20,1.10,1.02. 
The results in Fig. [^concern the sign of the second virial coefficient B 2 , red areas show 
a negative B 2 , where attractions overcome repulsions, and blue areas refer to a positive 
B 2 , where the repulsions are more important than the attractions. We remind the reader 
that the Boyle temperature kBT^oyXe/s is defined as the temperature where the second virial 
coefficient B 2 vanishes. We also note that the case of MM geometry, which matches the 
experimental geometry in Ref. [27], is found at the right-most end of each diagram, where 
l/ah = 0.57. 

For l/ah = 0.57, the region where R 2 < 0 is seen to dramatically shrink as the interaction 
range is decreased. This also means that the Boyle temperature is considerably dropping. 
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FIG. 7. The sign of the reduced second virial coefficient B 2 plotted in the reduced tem¬ 
perature k-QT/e versus bond length I/ah representation, for various interaction ranges A = 
1.50,1.40,1.30,1.20,1.10,1.02 as labelled. The orange triangle represent the state point where 
we have combined SUS and MC-NVT simulations as shown in Fig. |^and 10, the cyan triangles 
represent the state points we have investigated only via MC-NVT simulations, see Fig. [TOl 


going from k^T^oyie/s{\ — 1.50) ~ 0.9 to /cB2Boyie/£(A = 1.02) ~ 0.15. It is also known that 
temperatures below ~ 0.1 can get the system stuck in kinetic traps and our calculations show 
that upon decreasing the interaction range A up to A = 1.02 the Boyle temperature gets 
closer to this limit. In such conditions, the temperature window where GL phase separation 
can be found is very small, if present at all. A confirmation of our calculation is given in 
Ref. [27] , where it is mentioned that the temperature window between clustering and kinetic 
trapping is very narrow. 

Based on the second virial coefficient analysis, the largest temperature window where it 
is in principle still possible to locate a GL phase separation is A = 1.40, besides the already 
investigated case of A = 1.50. We have therefore calculated P{N) for a temperature range 
of k^T/e G [0.31,0.35], the resulting distributions being shown in Fig. 

For all investigated temperatures, the distributions appear to develop a shoulder at higher 
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particle number N than the the first peak corresponding to the gas phase. Upon increasing 
the system volume, this second peak is shown to remain at the same location, as shown 
in the inset of Fig. and is therefore associated with clusters developing in the system. 
This suggests that a self-assembly process preempts the macroscopic GL phase separation, 
either by shifting the transition to much lower temperatures than the investigated ones, 
or by destabilising the liquid with respect to the micellar fluid, as pointed out in previous 
works [HRl EDHHS]- Thus, even for the best case scenario of A = 1.40, we have not found 
evidence of a GL phase separation, but instead we find indication of spontaneously formed 
aggregates in the system. 



N 


- kj^T/e =0.35 - k^T/e =0.33 - k^T/e=0.3\ 


FIG. 8. Probability of finding N particles in the system for bond length I/ah = 0.57, interaction 
range A = 1.40 and temperature as labelled. All the distributions are reweighted as to highlight 
the peaked region present for particle number N G [10,40]. Inset: finite size study performed at 
temperature k^Tje = 0.31 confirming that the second peak is always located in the same window 
of particle number N. 


The snapshots of the system for particle numbers within the micellar peak of P(A^) 
(N = 15, N = 20, N = 20), shown in Fig. reveal the disordered microscopic configurations 
of the clusters associated with the relatively high interaction range A = 1.40. 

For lower values of the interaction range A, we select few state points, as indicated in Fig. 
and investigate the system via MG simulations in the NVT ensemble, due to computational 
reasons. The final configurations of the NVT simulations for A = 1.40,1.30,1.20,1.10 are 


reported in Fig. 10 The panel for a long ranged interaction A = 1.40 (top left) displays 
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FIG. 9. Disordered clusters of MMSW particles for bond length l/uh — 0.57, interaction range 
A = 1.40 and temperature k^T/e = 0.31. (left) N — 15, ((> = 0.011, (right) N — 20, cj) = 0.015. 

a micellar fluid with clusters of size comparable to the peak in the computed P{N). On 
the other hand, the remaining panels for shorter-ranged interactions show self-assembled, 
tube-like structures. Once the self-assembly regime has been reached, starting at A = 1.40, 
the shape of the resulting structures can be tuned by decreasing the interaction range. The 
backbone of the tubes becomes more close packed upon decreasing the interaction range, 
and the branching is suppressed for A < 1.05. Thus, structures as the Bernal spirals, which 
were reported in Ref. [27], are eventually recovered upon shrinking the interaction range to 
values lower than A = 1.05. 

The above findings highlight the role of the range in the self-assembly of MM particles: 
shrinking the interaction range first destabilises the GL phase separation with respect to 
self-assembly, and then favours the formation of increasingly ordered structures, eventually 
recovering those found in Ref. m- 

In a very recent work, the fluid phase behaviour of trimers with one soft, attractive bead 
and two soft, repulsive beads has been investigated [63] . There are substantial differences 
with our investigation, which we summarise in the following. Firstly, the particle model 
employed here uses hard-core particles, whereas the model in Ref. [63| makes use of soft¬ 
core particles. Secondly, the shape of their trimer is slightly different from our MM particle, 
since we keep the angle between the repulsive beads fixed. Thirdly, their methodology is 
fundamentally different, although it aims to compute the same quantity P{N). Lastly, due 
to our SW potential, we can explicitly control the interaction range between two particles, 
and track down the evolution of the self-assembly with it. Notwithstanding these differences, 
it is interesting to compare our results to Ref. |63| . Concerning the shift of the binodal line 
with the bond length distance, both works have found a linear scaling. The slopes of the 
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FIG. 10. Snapshots from MC NVT simulations at fixed packing fraction 0 = Nvp/V = 0.006 and 
bond length Ija^ = 0.57 for different interaction ranges A and temperatures k^T/e as described 
by the colored triangles in Fig. Top left: amorphous clusters for A = 1.40 and k^T/e — 0.310, 
in agreement with the P{N) calculated at the same temperature in Fig. Top right: branched 
tube-like structures appearing at A = 1.30 and k^T/e = 0.200; Bottom left: self-assembled tubes 
at A = 1.20 and k^Tje — 0.200; Bottom right: branched tubes at A = 1.10 and k^Tje — 0.143. 
For all the panels, N — 256 and V — 35( j ^. 

scaling are comparable, our estimate being a = —1.66, see the fit in Fig. against the 
reported a = —1.40. However, the actual values of the critical points for specific bond 
lengths do not match with each other, the discrepancy being ascribed to the differences 
in the particle geometry and particle-particle interaction. Furthermore, by very different 
analyses, both works agree on the region where the tube-like structures are expected to 
occur for the experiments reported in Ref. 123 - roughly for temperatures k-QT/e < 0.15. 


IV. CONCLUSIONS 

In this work we have investigated, via Monte Carlo simulations, how the combined ac¬ 
tions of particle geometry and particle-particle interaction range affect the gas liquid phase 
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separation. We have also addressed the role of the interaction range onto the self-assembly 
process. 

Starting from the well-known HSSW model, we have mapped out the different binodal 
curves for a family of Mickey Mouse particles with increasing steric repulsion. We have 
found that there is a linear relation between the critical temperature k-^Tc/s of the binodal 
curves and the bond length I/ah, which controls the influence of the repulsive spheres on 
the overall interaction. We interpreted this phenomenon in the “limited valence” paradigm 
and documented it by calculating the value of the second virial coefficient at the critical 
temperature B^^ik^Tc/e) as function of bond length l/au- This is a possible measure of the 
effective vafence of the system and indeed we have proven that this value decreases as the 
valence is reduced, which is controlled in this case via the bond length. 

The transition from macroscopic GL phase separation and self-assembly is driven by the 
interaction range, which has also an important effect on the Boyfe temperature Tsoyie of the 
system. In fact, as the interaction range is reduced, both the Boyle temperature decreases 
and the phase separation process is preempted by a self-assembly process. Microscopic con¬ 
figurations from the SUS simulations, complemented with snapshots from NVT simulations, 
have shown that the structure of the self assembled clusters changes from disordered to more 
ordered, branched tube-like structures. Upon further decrease in the range, these structures 
are seen to transform into the very ordered ones as reported in Ref. m- 

Our results highlight the role of particle geometry and interparticle interaction range 
in the competition between phase separation and self-assembly for Mickey Mouse particles 
and, in general, trimers with one attractive sphere and two repulsive spheres. This is of 
direct interest for experiments where trimers are used as basic building blocks, where - 
as we have documented - both the trimer exact shape and the trimer-trimer interaction 
range are expected to play a significant role. At a more general level, we have shown how 
tuning a particle’s valence can be achieved by controlling the bond length. This conclusion 
is expected to hold also in other particle geometries as long as the particles are composed of 
different beads, e.g. cofloidal dimers with one rough (repulsive) and one smooth (attractive) 
beads. 
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Appendix A: MM particle volume calculation 

The volume of a MM particle Vp can be calculated exactly via MC integration, and in an 
approximated fashion considering the amount of overlap volume between the two ears and 
the head. In fact, from geometrical considerations, it holds that 

Vp = Vh + 214 - 214,e - K,e (Al) 

where 14 is the volume of the attractive sphere, 14 the volume of the repulsive sphere, 
14,e is the head-ear overlap volume, and 14,e is the ear-ear overlap volume. Since, as first 
approximation, the two ears in a MM particle do not overlap with each other in the region 
outside the head, their overlap volume 14 ,e is small compared to the head-ear one 14 ,e so the 
latter is the only quantity we need to calculate, which is given by 

Vh,e = V{Rh,Re,l) = 

{Rh + Re “ 0 

= —^ 

12 / 

X + 2/Re - 3R^ + 2lRh + QReRh - ^Rl) (A2) 

where R^ = O.Sda, a = h, ei, 62 are the radii of the attractive and the repulsive spheres and 
/ is the center to center distance of the spheres. 

The Monte Carlo integration is similar to a MC calculation of tt, and we discuss it in the 
following. We place a MM particle in a box with fixed volume V, big enough to enclose the 
particle. We then generate a random position in the box and check whether it overlaps with 
the MM particle. If the number of generated random positions is Ntriai r\j (9(10®) and the 
number of hits on the MM particles is Nhu, then the MM particle volume is given by 

V, ~ (A3) 

trial 
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The results of the calculations for the volume of a MMSW particle are shown in Fig. [T^ 
for the ratio studied in this paper q = 0.85, as a function of bond length l/ah- The MC 
integration and the theoretical calculation are almost indistinguishable. 



FIG. 11. Volume of a MMSW particle as a function of bond length l/cih for size ratio q = 0.85. 
The solid green line is the prediction from Eq. |A1| with q — 0.85, and the blue dots are the results 


from the MC calculations according to Eq. |A3[ The errors are multiplied by 10 for visualisation 
purposes. 
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